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Abstract 



We study the time evolution of the entangled kaon system by considering 
the Liouville - von Neumann equation with an additional term which allows 
for decoherence. We choose as generators of decoherence the projectors to 
the 2-particle eigenstates of the Hamiltonian. Then we compare this model 
with the data of the CPLEAR experiment and find in this way an upper 
bound on the strength A of the decoherence. We also relate A to an effective 
decoherence parameter C, considered previously in literature. Finally we 
discuss our model in the light of different measures of entanglement, i.e. 
the von Neumann entropy S, the entanglement of formation E and the 
concurrence C, and we relate the decoherence parameter Q to the loss of 
entanglement: 1 — E. 
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I. INTRODUCTION 



Particle physics has become an interesting testing ground for fundamental questions of 
quantum mechanics (QM). For instance, QM versus local realistic theories [1-4] and Bell 
inequalities [5-12] have been tested. Furthermore, possible deviations from the quantum 
mechanical time evolution have been studied, particularly in the neutral K-meson system 
[13-19] and B-meson system [20-24]. Recently also neutrino oscillations have become of 
interest in this connection [25]. 

In this paper we concentrate on possible decoherence effects arising due to some 
interaction of the system with its "environment". Sources for "standard" decoherence 
effects are the strong interaction scatterings of kaons with nucleons, the weak interaction 
decays and the noise of the experimental setup. "Nonstandard" decoherence effects result 
from a fundamental modification of QM and can be traced back to the influence of 
quantum gravity [26-28] - quantum fluctuations in the space-time structure on the Planck 
mass scale - or to dynamical state reduction theories [29-32], and arise on a different 
energy scale. However, we do not pursue further the reasons for decoherence effects, rather 
we want to develop a specific model of decoherence and quantify the strength A of such 
possible effects with the help of data of existing experiments. 

For our model we focus on entangled massive particles moving apart in their center 
of mass system, in particular on the K^K^ system, where the strangeness S = 
plays the role of spin "up" and "down" (for details see Ref. [33]). We consider here the 
famous EPR-like (Einstein, Podolsky. Rosen) scenario, as described by Bell [34] for spin- 
1/2 particles, where the initial spin singlet state evolves in time and after macroscopic 
separation the strangeness of the left and right moving particle is measured. In contrast 
to other concepts in the literature we introduce decoherence in the time evolution of the 
2-particle entangled state which becomes stronger with increasing distance between the 
two particles, whereas for the 1-particle state we assume the usual quantum mechanical 
time evolution. 

Then we compare our model of decoherence with the experimental data of the 
CPLEAR experiment performed at CERN [35] and find an upper bound on possible 
decoherence. We also can relate our model to an effective decoherence parameter ^, in- 
troduced previously in literature, which quantifies the spontaneous factorization of the 
wavefunction into product states (Furry-Schrodinger hypothesis [36,46]). 

Finally we discuss our model within concepts of quantum information, where en- 
tanglement is quantified by certain measures. We can connect directly the amount of 
decoherence of the K^K^ system parametrized by A or with the loss of entanglement 
expressed in terms of the concurrence or in terms of entanglement of formation. The 
numerics for the information loss, the von Neumann entropy, and the entanglement loss 
of the evolving K°K° system we illustrate in Fig.l. 
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II. THE MODEL 



Let us begin our decoherence discussion with the l-particle kaon system as an intro- 
duction. Then we proceed to the case of two entangled neutral kaons and compare it with 
experimental data. 

A. The l-particle case 

We discuss the model of decoherence in a 2-dimensional Hilbert space Ti = and 
consider the usual non-Hermitian "effective mass" Hamiltonian H which describes the 
decay properties and the strangeness oscillations of the kaons. The mass eigenstates, the 
short hved \Ks) and long lived \Ki) states, are determined by 

H\Ks,l) = >^s,l\Ks,l) with Xs,l = rns,L - s,l , (2.1) 

with ms,L and Ts,l being the corresponding masses and decay widths. For our purpose 
CV invariance^ is assumed, i.e. the CV eigenstates \Ki), \K2) are equal to the mass 
eigenstates 

\Kl) = \Ks), \KI) = \Kl), and {Ks\Kl) = Q. (2.2) 

As a starting point for our model of decoherence we consider the Liouvillc - von 
Neumann equation with the Hamiltonian (2.1) and allow for decoherence by adding a 
so-called dissipator D[p], so that the time evolution of the density matrix p is governed 
by a master equation of the form 

^ = -iHp + ipH^ -D[p]. (2.3) 
For the term D[p] we choose the following ansatz (as in Ref. [37]) 

D[p] = \{PspPl + PlpPs) = \Y1 \P3^[P3.P]] , (2-4) 

3=S,L 

where Pj = \Kj){Kj\ {j = S,L) represent the projectors to the eigenstates of the 
Hamiltonian and the decoherence parameter A is positive, A > 0. 

Apart from its simplicity ansatz (2.4) has the following nice features: 



^Note that corrections due to CV violations are of order 10~^, however, we compare this 
model of decoherence with the data of the CPLEAR experiment [35] which are not sensitive to 
CV violating effects. 
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i) It generates a completely positive map since it is a special case of Lindblad's general 
structure [38] 

D\P\ = \ Y.^A]A, p + pA]A, - 2A,pA]) (2.5) 
j 

if we identify Aj = \f\Pj , j = S, L. 

Equivalcntly, it is a special form of the Gorini-Kossakowski-Sudarshan [39] expres- 
sion (see, e.g., Ref. [40]). 

ii) It conserves energy in case of a Hermitian Hamiltonian since [Pj,H] — (see, e.g. 
Ref. [41]). 

iii) The von Neumann entropy S{p) — — Tr(plnp) is not decreasing as a function of 
time since P/ = Pj, thus A] = A , in our case, which is a theorem due to Narnhofer 
and Benatti [42]. 

With choice (2.4) the time evolution (2.3) decouples for the components of p which 
are defined by 

P(t) = Yl P^At)\K^){K,\, (2.6) 
i,j=S,L 

and we obtain 

Pssit) = pssiO) ■ e-^^* 
PLL(t) = PLL(0)-e-r^* 

PLsit) = PLsiO) ■ e-^-*-r*-^* (2.7) 

with Am = rriL — mg and F = \{ys + T/,). Only the off-diagonal elements are effected 
by our model of decoherence. Before discussing further the model we now proceed to the 
2-particle system. 

B. The 2-particle case 

In the case of two entangled neutral kaons we make the following identification 

|ei) = \Ks)i®\KL)r and jes) = \KL)i®\Ks)r , (2.8) 

and we consider - as common - the total Hamiltonian as a tensor product of the 1- 
particle Hilbert spaces: if = if; ® 1^ + 1; ® if^, where / denotes the left-moving and r 
the right-moving particle. Then the initial quasispin singlet state 

lV^-) = ^{|ei)-|e2)}, (2.9) 
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is equivalently given by the density matrix 



p(0) = = -{|ei)(ei| + |e2)(e2|-|ei)(e2|-|e2)(ei|}. (2.10) 

Following the considerations of the 1-particle case wc find that the time evolution given 
by (2.3) with our choice (2.4), where now the operators Pj = \ej){ej\ {j — 1,2) project 
to the eigenstates of the 2-particle Hamiltonian H, also decouples 

Pii(t) =Pii(0) e-'^* 
P22(t) =P22(0) e-2r* 

Pi2(i)=Pi2(0)e-2r*-^*. (2.11) 
Consequently we obtain for the time-dependent density matrix 

Pit) = ^e-2r*||ei)(ei| + |e2)(e2| - e-^*( |ei)(e2| + |e2)(ei| )} . (2.12) 

The decoherence arises through the factor e~'^* which only effects the off-diagonal ele- 
ments. It means that for t > and A 7^ the density matrix p{t) does not correspond to 
a pure state anymore (for further discussion, see Section IV). 

Note that the assumption of CV invariance, Eq.(2.2), - which is sufficient for our 
purpose - implying (ei|e2) = 0, is crucial. Otherwise we would have a time evolution into 
the full 4-dimensional Hilbert space of states. 



C. Bounds from experimental data 

In order to obtain information on possible values of A we compare our model of 
decoherence with data of the CPLEAR experiment performed at CERN [35]. We have 
the following point of view. The 2-particle density matrix follows the time evolution given 
by Eq.(2.3) with the Lindblad generators Aj = \/X |ej)(ej| and undergoes thereby some 
decoherence. We measure the strangeness content S of the left-moving particle at time 
ti and of the right- moving particle at time tr- For sake of definiteness we choose U < ti. 
For times tr < t < ti we have a 1-particle state which evolves exactly according to QM, 
i.e. no further decoherence is picked up. 

Mathematically, the measurement of the strangeness content, i.e. the right-moving 
particle being a or a at time tr, is obtained by 

Tr,{l, \S){S'\r P{tr)} = Pl{t = tr\tr) , (2.13) 

with strangeness S' — -\-,— and |-|-) = \K^), |— ) = 1-^°)- Consequently, pi{t;tr) for times 
t > tr is the 1-particle density matrix for the left- moving particle and evolves as a 1- 
particle state according to pure QM. At t = ti the strangeness content (5* = -|-, — ) of the 
second particle is measured and we finally have the probability 

P^{S,ti;S\tr)^T:ii{\S){S\i pi{ti;tr)} . (2.14) 
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Explicitly, we find the following results for the like- and unlike-strangeness probabilities 
Px{K^,ti\K'',tr) = Px{K',ti;K'',tr) = 

= 1 |e-rs*.-r^t. + g-r^t.-r^t. _ ^-xu 2cos(AmAt) • e-^(*'+*^)} 

= ^ je-'^s*'-''^*'- + e-^^*'-''s*r ^ g-At. 2 cos( AmAi) • e-''(*'+*^) | , (2.15) 

with At — ti — tr- 

Note that at equal times ti—tr — t the like-strangeness probabihties 

Px{K',t;K',t) = Px{K\t;K\t) = 1 e'^^* (1 - e"^*) (2.16) 
do not vanish, in contrast to the quantum mechanical EPR-correlation. 

The asymmetry of the probabilities is directly sensitive to the interference term and 
has been measured by the CPLEAR collaboration [35]. For pure QM we have 

AQ^^At) = 

P{K^, tu K^, tr) + P{K^, tu K^, tr) - PiK^, tu K^, tr) - P{K\ tu K^, tr) 



P{KO, tr, K^, tr) + P{KO, ti- K^, tr) + P{K^, ti- XO, tr) + P{K\ t^ K\ tr) 

cos(AmAt) 
cosh(iArAt) ' 



(2.17) 



with Ar = Tl — Tsi and for our decoherence model we find by inserting the probabilities 
(2.15) 

A'^itutr) = ^Qg^tr^l) ■ e-^""^*"*"^^ = A^^{At) ■ e-^-'° . (2.18) 

Thus the decoherence effect, simply given by the factor Q-X'^^itiM ^ depends only - due 
to our philosophy - on the time of the first measured kaon, in our case: min {ti, tr} — tr- 

Comparing now our model with the results of the CPLEAR experiment [35] we recall 
that the experimental set-up has two configurations: In the first configuration both kaons 
propagate 2cm, in the second one kaon propagates 2cm and the other kaon 7cm until 
they are measured by an absorber. 

Fitting the decoherence parameter A by comparing the asymmetry (2.18) with the 
experimental data [35] we find, when averaging over both configurations^, the following 
bounds on A 



^We have scaled At in the QM asymmetry (2.17) in order to reproduce the QM curve in Fig. 
9 of the CPLEAR collaboration [35]. 
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A = (1.84l|f°) • 10-^2 MeV and A = — = 0.25l°|^ . (2.19) 

The results (2.19) are certainly compatible with QM (A = 0), nevertheless, the experi- 
mental data allow an upper bound Xup — 4.34 • 10~^^ MeV for possible decoherence in 
the entangled K^K^ system. 

The results (2.19) can be compared with the analogous investigation of the entangled 
system [37], which gives Xb = (-47 ± 76) • lO'^^ ^eV. Thus we find that the 
bounds (2.19) of the K^K^ system are an order of magnitude more restrictive. 

III. CONNECTION TO A PHENOMENOLOGICAL MODEL 

There exists a one-to-one correspondence between the model of decoherence (2.3) and 
a phenomenological model [22,43-45] where the decoherence is introduced by multiplying 
the interference term of the transition amplitude by the factor (1 — Q. The quantity C, 
is called the effective decoherence parameter. Our initial state is again the spin singlet 
state which is given by the mass eigenstate representation (2.9) then we get for the 
like-strangeness probability 

PiK',ti;K',tr) = \\{K\^{KW^Ij-iti,tr))\\'' PciK',ti;K'',tr)^ 
modification 

= l|e-rst<-r.t. ^ g-r,t,-r,t. _ 2 (^-C) cos(AmAO • e-r(*'+*^)} , 

modification (3.1) 

and the unlike-strangeness probability just changes the sign of the interference term. 

The value ( = corresponds to pure QM and C = 1 to total decoherence or sponta- 
neous factorization of the wave function, which is commonly known as Furry-Schrodinger 
hypothesis [36,46]. The effective decoherence parameter introduced in this way "by 
hand" , interpolates continuously between these two limits and represents a measure for 
the amount of decoherence which results in a loss of entanglement of the total quantum 
state (we come back to this point in Section IV). 

Calculating the asymmetry of the strangeness events with the probabihties (3.1) we 
obtain 

A'^{ti,tr) = A'^^\At)-{l-ati,tr)). (3.2) 

When we compare now the two approaches, i.e. Eq.(2.18) with Eq.(3.2), we find the 
formula 
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(3.3) 



Of course, the values (2.19) are in agreement with the corresponding ( values (averaged 
over both experimental setups): ( — 0.131q;}5 , as derived in Refs. [43,47]. 

We consider the decoherence parameter A to be the fundamental constant, whereas the 
value of the effective decoherence parameter ( depends on the time when a measurement 
is performed. In the time evolution of the state IV'"), Eq.(2.9), represented by the density 
matrix (2.12), we have the relation 

at) = 1 - e-'' , (3.4) 

which after the measurement of the left- and right moving particles at ti and U turns 
into formula (3.3), when decoherence is implemented at the 2-particle level. 

Our model can be compared with the case where decoherence in the time evolution 
(2.3) happens at a 1-particle level and is transferred to the 2-particle level by a tensor 
product of the 1-particle Hilbert spaces [18]. Using the same structure of the decoherence 
term (2.4), where now the operators project to the 1-particle states instead of states (2.8), 
we obtain the relation 

ati,tr) = l-e-^(*'+*'^). (3.5) 

Here the parameter ( depends explicitly on both times ti and tr, the "eigentimes" of the 
left- and right-moving kaon, instead of one time min {ti, tr}, the time of the system when 
the first kaon is measured. 

Measuring the strangeness content of the entangled kaons at definite times we have 
the possibility to distinguish experimentally these two models (3.3) and (3.5) on the basis 
of time dependent event rates. Indeed, it would be of high interest to measure in future 
experiments the asymmetry of the strangeness events for several different times, in order 
to confirm the time dependence of the decoherence effect. In fact, such a possibility is 
now offered in the B-meson system. Entangled pairs are created with high density 

at the asymmetric B-factories and identified by the detectors BELLE at KEK-B (see e.g. 
Refs. [48,49]) and BABAR at PEP-II (see e.g. Refs. [50,51]) with a high resolution at 
different distances or times. 

IV. DECOHERENCE AND ENTANGLEMENT LOSS 

The term D[p] in the master equation (2.3) is usually called dissipative term or dissi- 
pator (see, e.g., Ref. [52]). In general, D[p] describes 2 phenomena occurring in an open 
quantum system S, namely decoherence and dissipation. When the system S interacts 
with the environment E the initially product state evolves into entangled states of S + E 
in the course of time [53,54]. It leads to mixed states in S - what means decoherence - 
and to an energy exchange between S and E - what is called dissipation. 
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The decoherence destroys the occurrence of longrange quantum correlations by sup- 
pressing the off-diagonal elements of the density matrix in a given basis and leads to an 
information transfer from S to E. 

In general, both effects are present, however, decoherence acts on a much shorter time 
scale [53,55-57] than dissipation and is the more important effect in quantum information 
processes. 

Our model describes decoherence and not dissipation. The increase of decoherence of 
the initially totally entangled K^K^ system as time evolves means on the other hand a 
decrease of entanglement of the system. This loss of entanglement we can quantify and 
visualize explicitly. 

In the field of quantum information the entanglement of a state is quantified by 
introducing certain measures. In this connection the entropy plays a fundamental role, 
which measures "somehow" the degree of uncertainty of a quantum state. A common 
measure is the von Neumann entropy function S, the entanglement of formation E and 
the concurrence C. 



A. Von Neumann entropy 

Since we arc only interested in the effect of decoherence, we want to properly nor- 
malize the state (2.12) in order to compensate the decay property^ of the non-Hermitian 
Hamiltonian H 

Mt) = ^ . (4.1) 

Then von Neumann's entropy function for the state (4.1) gives 

S{pM{t)) = -Tr{pr,{t)\og^pN{t)} 

l-e-^\ 1-e-^* 1 + 1 + e"^* , , 

= ^ log2 2— log. . (4.2) 

At the time t — the entropy is zero, there is no uncertainty in the system, the quantum 
state is pure and maximally entangled. For t > the entropy gets nonzero, increases 
and approaches the value 1 for t — > oo. Hence the state becomes more and more mixed. 
Mixed states provide only partial information about the system, and the entropy mea- 
sures how much of the maximal information is missing (see, e.g., Ref. [58]). In Fig.l the 
von Neumann entropy S'(pjv(t)) is plotted for the mean value and upper bound of the 
decoherence parameter A, Eq.(2.19), as determined from the CPLEAR experiment [35]. 

Let us consider next the reduced density matrices of the subsystems, i.e. the propa- 
gating kaons on the left I and right r hand side 



^Note that for other physical situations - like the verification of Bell inequalities - the decay 
property must not be neglected [8] . 
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pj,{t)^T:tr{pN{t)} and p'^{t)^T:ti{pN{t)} 



(4.3) 



Then the uncertainty in the subsystem I before the subsystem r is measured is given by 
the von Neumann entropy S{pj^{t)) of the corresponding reduced density matrix plf{t) 
(and alternatively we can replace I ^ r). In our case we find 

sipUt))^sip^^it))^i yt>o. (4.4) 

We see that the reduced entropies are independent of A. The correlation stored in the 
composite system is, with increasing time, lost into the environment - what is expected 
intuitively - and not into the subsystems, i.e. the individual kaons. 

Note that because of our chosen normalization (4.1) the effects seen arc only due 
to the introduced decoherence via the term D[p] (2.4) and not due to the decay of the 
system. 

B. Lack of separability 

For the subsequent considerations it is convenient to recall the "quasispin" picture for 
the K^K'^ system (see, e.g., Ref. [8]) in order to express the formulae in terms of Pauli 
spin matrices. Then the projection operators to the mass eigenstates represent the spin 
projection operators "up" and "down" 

Ps^\Ks){Ks 

and the transition operators are the "spin-ladder" operators 

Psl = \Ks){Kl 
Pls^\Kl){Ks 

With help of these spin matrices the density matrix (4.1) can be expressed by 

PN{t) = i{(T| ^a^ + a^^a^ - e"^* [a+ ^ a_ + a_ ^ a+]} , (4.7) 

or by 

PAf(^) = ^{1 - (7z®(yz - e"^* {Ox (7a; + (7j, CTj/]} , (4.8) 

which at f = coincides with the well-known expression for the pure spin singlet state 
pjsjit = 0) = I {\ — a®a)\ see, e.g., Ref. [59]. Operators (4.7), (4.8) can be nicely written 
as 4 X 4 matrix 



1 





1 



(4.5) 



0-+ = -^{(^x + iCTy) 



-[a^-t ay) 



1 




1 



(4.6) 
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PN{t) 





/o 








o\ 


1 





1 







2 







1 







^0 








0/ 



(4.9) 



It is also illustrative to choose for the density matrix pjv(t) another basis, the socalled 
"Bell basis" 



— 1-0^) (-0^1 and 
with 1-0") given by Eq.(2.9) and 1-0+) by 

l^->-^{h) + |e.>}. 



(4.10) 



(4.11) 



The states = -^(| tt) =F I ii)) (in spin notation) do not contribute here. Then we 
are led to the following proposition. 

Proposition: The state represented by the density matrix pN{t) (4-1) becomes mixed for 

< t < cxD hut remains entangled. Separability is achieved asymptotically t oc with the 
weight e~^^ . Explicitly, pN{t) is the following mixture of the Bell states p~ and p"*" : 



PN{t) 



-At 



)P- 



(4.12) 



Proof: 



i) Diagonalizing matrix (4.9) we find + e""^*), |(1 — e""^*), 0, 0}, the eigenvalues 
together with the eigenvectors and p+ , which proves Eq.(4.12). For t = the 
matrix pAr(O) is a projector to a 1-dimcnsional subspace: {1,0,0,0}, the Bell state 
p~ , so the state is pure. For < t < cxo the matrix pAr(t) is no longer a projector, 
thus the state is mixed. The state becomes asymptotically t oo totally mixed, 
separable. Of course, pN{t) also satisfies the mixed state criterion 



put) 





/o 








o\ 


1 





1 + 


-2e-^* 





4 





-2e-^* 


1 + e-2At 







^0 








0/ 



^ PN{t) for i > 



(4.13) 



ii) Entanglement we prove via lack of separability. According to Peres [60] and the 
Horodeckis [61] separability is determined by the positive partial transposition cri- 
terion: the partial transposition of a separable state with respect to any subsystem 
is positive (i.e. the operator has positive eigenvalues). Applying the transposition 
operator T, defined by T{a'^)ki — {a'^)ik , to the left- or right hand side, we find a 
negative eigenvalue: 
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(1; (g) Tr) pN{t) = ^{u^ ^ai + ai^a^ - e ^* [a+ ® (7+ + (7_ ® (7_]} 



/ -e-^*\ 




\ -e"^* / 



1 

1 
e"^* 



^0, 



(4.14) 



2 



with eigenvalues {2,2)2^ 

Alternatively, we could use Horodecki's [62] reduction criterion: pN{t) is separable 
iff Pj^it) ® — pN{t) > or 1; ® p^(t) — pN(t) > 0. In our case we have 



h ^ PN{t) - PNit) 



/I 0\ 

e-^* 
e-^* 

yo 01/ 



^0, 



(4.15) 



where the eigenvalues are the same as for the previous criterion. For a general 
separability criterion, a generalized Bell inequality, see Ref. [59] . □ 



C. Entanglement of formation and concurrence 

For pure quantum states entanglement can be measured by the entropy of the reduced 
density matrices. For mixed states, however, the von Neumann entropy is not generally 
a good measure for entanglement, whereas another measure, entanglement of formation 
[63], is very suitable. It has been constructed by the authors [63] to quantify the resources 
needed to create a given entangled state. 



1. Definitions 

Every density matrix p can be decomposed in an ensemble of pure states pi — 
with the probability pi, i.e. p = J2iPiPi- The entanglement of formation for a pure state is 

given by the entropy of either of the two subsystems. For a mixed state the entanglement 
of formation for a bipartite system is then defined as the average entanglement of the 
pure states of the decomposition, minimized over all decompositions of p 

E{p) ^ mm J2 Pi S{p^). (4.16) 

i 

Bennett et al. [63] found a remarkable simple formula for entanglement of formation 

E{p) > 8{f{p)) , (4.17) 
where the function S{f{p)) is defined by 
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mP)) - + v/7(W)) for / > i , (4.18) 

and S{f{p)) = for / < |. The function H represents the familiar binary entropy function 
H{x) = — xloga^ — (1 — x) log2(l — a;). The quantity / is called the fully entangled fraction 
of p 

f = max(e|p|e) , (4.19) 

being the maximum over all completely entangled states |e). For general mixed states 
p the function £(f(p)) is only a lower bound to the entropy E{p). For pure states and 
mixtures of Bell states - the case of our model - the bound is saturated, E — £, and we 
have formula (4.18) for calculating the entanglement of formation. 

Wootters and Hill [64-66] found that entanglement of formation for a general mixed 
state p of two qubits can be expressed by another quantity, the concurrence C 

E{p)^S{C{p)) = + wi^h 0<C<1. (4.20) 

Explicitly, the function £{C) looks like 



^ iC) = 7\ log2 7, 7, log2 7 



(4.21) 



and is monotonically increasing from to 1 as C runs from to 1. Thus C itself is a kind 
of entanglement measure in its own right. 
Defining the spin flipped state p of p by 

p= {(Ty® ay) p*{ay (g) ay) , (4.22) 

where p* is the complex conjugate and is taken in the standard basis, i.e. the basis 
{| TT)) I ii)) I Ti)) I it)}) the concurrence C is given by the formula 

C{p) = max{0, Ai - A2 - A3 - A4} . (4.23) 

The Aj's are the square roots of the eigenvalues, in decreasing order, of the matrix pp. 

2. Applications to our model 

For the density matrix pAr(^) (4.1) of our model, which is invariant under spin flip 
(see, e.g., Eq.(4.7)), i.e. pjv = Pn and thus PnPn — Pn, we get for the concurrence 

C{pN{t)) =max{0,e-^'} = e'^\ (4.24) 

and for the fully entangled fraction of pN{t) 
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/(PivW) = ^(l + e-^*), (4.25) 

which is simply the largest eigenvalue of pNif)- Clearly, in our case the functions C and 
/ are related by 

C{pM{t))^2f{pM{t))-l. (4.26) 
Finally, for the entanglement of formation of the K^K^ system we have 



E{pN{t)) = log2 



e 



-2 At 



1 - Vl - e-2At 1- 
2 ^""^^ 2 ■ ^^-^^^ 

Recalling our relation (3.4) between the decoherence parameters A and C, we find a di- 
rect connection between the entanglement measure and the amount of decoherence of 
the quantum system. Defining the loss of entanglement as one minus entanglement and 
expanding expression (4.27) for small values of A or C we obtain 

1 - C{pN{t)) = C{t) , (4.28) 
1 - E{p^{t)) = ^ C(^) = 1^ ^ ■ (4-29) 

The loss of entanglement of the propagating K^K° system in terms of the concurrence 
C, Eq.(4.28), equals precisely the amount of the decoherence parameter ( which describes 
the factorization of the initial spin singlet state into the product states \Ks)i ^ \KL)r or 
\Kl)i ® \Ks)r (Furry-Schrodingcr hypothesis). In terms of entanglement of formation, 
Eq.(4.29), the decoherence parameter is weighted by a factor l/ln2 = 1.44 (the factor 
In 2 reflects the dimension 2 of the X-meson quasispin space) and the entanglement loss 
increases linearly with time. 



D. Discussion of the results 



In Fig.l we have plotted the loss of entanglement 1 — E, given by Eq.(4.27), as 
compared to the loss of information, the von Neumann entropy function S, Eq.(4.2), in 
dependence of the time t/Tg of the propagating K^K^ system. The curves are shown for 
the mean value and upper bound of the decoherence parameter A, Eq.(2.19). The von 
Neumann entropy function visualizes the loss of the information about the correlation 
stored in the composite system. Remember that the information flux is not flowing into 
the subsystems but into the environment, see Eq.(4.4). The loss of entanglement of for- 
mation increases slower with time and visualizes the resources needed to create a given 
entangled state. At i = the pure Bell state p~ is created and becomes mixed for t > 
by the other Bell state p"*". In the total state the amount of entanglement decreases until 
separability is achieved (exponentially fast) for t — > oo. 
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FIG. 1. The time dependence of the von Neumann entropy (dashed hnes), Eq.(4.2), and 
the loss of entanglement of formation 1 — E (solid lines), given by Eq.(4.27), are plotted 
for the experimental mean value A = 1.84 • 10~^^ MeV (lower curve) and the upper bound 
Xup = 4.34 • 10-12 (upper curve), Eq.(2.19), of the decoherence parameter A. The time t 

is scaled versus the lifetime Tg of the short lived kaon Ks- t — > t/Ts- The vertical lines represent 
the propagation time Iq/ts ~ 0.55 of one kaon, including the experimental error bars, until it 
is measured by the absorber in the CPLEAR experiment. 

For example, in case of the CPLEAR experiment, where one kaon propagates about 
2 cm, which corresponds to a propagation time Iq/ts ^ 0.55, until it is measured by an 
absorber, the loss of entanglement is about 18% for the mean value and maximal 38% for 
the upper bound of the decoherence parameter A. These values, however, could diminish 
considerably in future experiments. 



We have considered a simple model of decoherence of the entangled K^K^ state due to 
some environment, i.e. a master equation of Liouville - von Neumann type with an addi- 
tional term D[p\. As generators of D[p\ causing the decoherence effect with strength A we 
choose the projectors to the eigenstates of the "effective mass" Hamiltonian and, for sim- 
plicity, CV invariance is assumed. For this choice the time evolution for the components 
of p decouples and only the off-diagonal elements are effected by our modification. 

We apply the model to the data of the CPLEAR experiment, where we follow the 
philosophy that only the 2-particle state is affected by decoherence, whereas the 1-particle 
state evolves according to pure QM. We estimate in this way the strength A of the 
occurring decoherence, Eq.(2.19). 

Moreover, we can relate the model to the case of a phenomenologically introduced 
decoherence parameter C, and find a one-to-one correspondence, Eq.(3.4). However, the 



V. SUMMARY AND CONCLUSIONS 
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existing data are not yet sufficient to measure the time-dependence of ( as predicted by 
our model, Eq.(3.3). So further measurements of the time-dependent asymmetry term, 
Eqs.(2.18), (3.2), would be of high interest in future experiments and will sharpen con- 
siderably the bounds of the parameters A and 

We can directly relate the decoherence of the K^K^ state to its loss of entanglement. 
In this connection we consider entanglement measures frequently discussed in the field 
of quantum information. We demonstrate that the initially pure singlet state of the 
entangled K^K^ becomes mixed for {) < t < oo but remains entangled and achieves 
separability for i — > oo (see Proposition). 

We find that entanglement loss in terms of the concurrence equals precisely the de- 
coherence parameter Eq.(4.28), and in terms of entanglement of formation the loss 
is very well approximated by or ■ t, Eq.(4.29), which is one of our main results. 
We can propose in this way how to measure experimentally the entanglement of the 
K^K^ system. 

In Fig.l we visualize both the loss of information given by the von Neumann entropy 
- which flows totally into the environment and not into the subsystems of the 2-particle 
system - and the loss of entanglement of formation. Inserting the mean value and upper 
bound of the parameter A, which we have determined from the CPLEAR experiment, we 
obtain definite bounds for both the information- and entanglement loss of the propagating 
K^K^ system in dependence of the time. These values, however, could be improved 
considerably in future experiments. 
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